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We study the flow induced by random vibration of a solid boundary in an otherwise 
quiescent fluid. The analysis is motivated by experiments conducted under the low level 
and random effective acceleration field that is typical of a microgravity environment. 
When the boundary is planar and is being vibrated along its own plane, the variance 
of the velocity field decays as a power law of distance away from the boundary. If a 
low frequency cut-off is introduced in the power spectrum of the boundary velocity, the 
variance decays exponentially for distances larger than a Stokes layer thickness based on 
the cut-ofF frequency. Vibration of a gently curved boundary results in steady streaming 
in the ensemble average of the tangential velocity. Its amplitude diverges logarithmically 
with distance away from the boundary, but asymptotes to a constant value instead if 
a low frequency cut-off is considered. This steady component of the velocity is shown 
to depend logarithmically on the cut-off frequency. Finally, we consider the case of a 
periodically modulated solid boundary that is being randomly vibrated. We find steady 
streaming in the ensemble average of the first order velocity, with flow extending up to a 
characteristic distance of the order of the boundary wavelength. The structure of the flow 
in the vicinity of the boundary depends strongly on the correlation time of the boundary 
velocity. 



1. Introduction 

This paper examines the formation and separation of viscous layers in a fluid which 
is in contact with a solid boundary that is vibrated randomly. The analysis is motivated 
by the low level and random acceleration field that affects a number of microgravity 
experiments. We first study the case of a planar boundary to generalize the classical 



result of Stokes (1851) who considered a boundary vibrated periodically along its own 
plane. We next consider a slightly curved boundary, and show that steady streaming 
appears in the ensemble average at first order in the perturbed flow variables. There 



are several qualitative similarities and differences with the classical result by Schlichting 



(1979) for the case of periodic vibration. Finally, we address the case of a modulated 
boundary that is vibrated randomly. 

Our study is motivated by the significant levels of residual accelerations (g-jitter) that 
have been detected during space missions in which microgravity experiments have been 
conducted ( [Walter (1987)| , Nelson (1991), DcLombard et al. (1997)). Direct measurement 



of these residual accelerations has shown that they have a wide frequency spectrum, 
ranging approximately from 10~'^Hz to lO^Hz. Amplitudes range from 10~^gE at the 
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lowest end of the frequency spectrum, and increase roughly linearly for high frequencies, 
reaching values of lO^^gg — \Q~^gE at frequencies of around IQHz {gs is the intensity 
of the gravitational field on the Earth's surface). Despite the efforts of a number of 
researchers over the last decade, there remain areas of uncertainty about the potential 
effect of such a residual acceleration field on typical microgravity fluid experiments, 
especially in quantitative terms. A better understanding of the response of a fluid to such 
disturbances would enable improved experiment design to minimize or compensate for 
their influence. In addition, it would also be useful to have error estimates of quantities 
measured in the presence of residual accelerations, including whenever possible some 
methodology for extrapolation to ideal zero gravity. 

The formation of viscous layers around solid boundaries when the flow amplitude 
has a random component has not been addressed yet despite its potential relevance 
for a number of microgravity experiments. Among them we mention the dynamics of 
colloidal suspensions, coarsening studies of solid-liquid mixtures in which purely diffusive 
controlled transport is desired, or the interaction between the viscous layer produced 
by bulk flow of random amplitude and the morphological instability of a crystal-melt 
interface. Our study represents the flrst step in this direction, and focuses on simple 
geometries in order to elucidate those salient features of the flow that arise from the 
random nature of the vibration. 

Previous theoretical work on the influence of g-jitter on fluid flow ranges from order 
of magnitude estimates to detailed numerical calculations. For example, the order of 
magnitude of the contributions to fluid flow from the residual acceleration field may be 
estimated from the length and time scales of a particular experiment, and the values of the 



relevant set of dimensionless numbers (Alexander (1990)). Such studies are of interest as 



a first approximation, but are not very accurate. Other studies have modeled the residual 
acceleration field by some simple analytic function in which the acceleration is typically 
decomposed into steady and time dependent components, the latter being periodic in 



time (Gershuni & Zhukhovitskii (1976), Kamotani et al. (1981), Alexander et al. (1991 
Farooq fc Homsy (1994) , Grassia & Homsy (1998a), prassia fc Homsy (1998&) , Gershuni 



fc Lyubimov (1998)1 ). A few studies have also a ddressed the consequences of isolated 
pulses of short duratio n ( [Alexander et al. (1997) ). 

/jtiang et al. (lyyaj and Thomson et al. (1997) adopted a statistical description of the 



residual acceleration field onboard spacecraft, and modeled the acceleration time series 
as a stochastic process in time. The main premise of this approach is that a statistical 
description is necessary in those cases in which the characteristic time scales of the phys- 
ical process under investigation are long compared with the correlation time of g-jitter, 
T (the acceleration amplitudes and orientations at two different times are statistically in- 
dependent if separated by an interval larger than r) . Progress has been achieved through 
the consideration of a specific stochastic model according to which each Cartesian com- 
ponent of the residual acceleration field g{t) is modeled as a narrow band noise. This 
noise is a Gaussian process defined by three independent parameters: its intensity (5^), 
a dominant angular frequency fi, and a characteristic spectral width r^^. Each realiza- 
tion of narrow band noise can be viewed as a temporal sequence of periodic functions 
of angular frequency with amplitude and phase that remain constant only for a fi- 
nite amount of time (r on average). At random intervals, new values of the amplitude 
and phase are drawn from prescribed distributions. This model is based on the following 
mechanism underlying the residual acceleration field: one particular natural frequency 
of vibration of the spacecraft structure (il) is excited by some mechanical disturbance 
inside the spacecraft, the excitation being of random amplitude and taking place at a 
sequence of unknown (and essentially random) instants of time. 
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Narrow band noise has been shown to describe reasonably weU many of the features 



of g-jitter time series measm^ed onboard Space Shuttle by Thomson et al. (1997). Actual 
g(-jitter data collected during the SL-J mission were analyzed, and a time series of roughly 
six hours was studied in detail. A scaling analysis revealed the existence of both deter- 
ministic and stochastic components in the time series. The deterministic contribution 

1/2 

appeared at a frequency of 17 Hz, with an amplitude (g^) — 3.56 x lO^^gE- Stochas- 
tic components included two well defined spectral features with a finite correlation time; 
one at 22 Hz with (g^)^^^ = 3.06 x IQ-^gs and r = 1.09 s, and one at 44 Hz with 

{g'^Y^^ — 5.20 X lO^'^gE and r = 0.91 s. White noise background is also present in the 
series with an intensity D — 8.61 x lO^^cm^/s"^. 

A further theoretical advantage of narrow band noise is that it provides a convenient 
way of interpolating between monochromatic noise (akin to studies involving a determin- 
istic and periodic gravitational field) , and white noise (in which no frequency component 
is preferred). In the limit r — s- with D = (^g"^^ r finite, narrow band noise reduces to 
white noise of intensity Z?; whereas, for r — )■ oo with (5^) finite, monochromatic noise is 
recovered. 

We discuss in this paper the flow induced in an otherwise quiescent fluid by the random 
vibration of a solid boundary. The velocity of the boundary Mo(i) is assumed prescribed, 
and modeled as a narrow band stochastic process. First, we consider an inflnite planar 
boundary that is being vibrated along its own plane to generalize the classical problem 



studied by Stokes (1851). In the monochromatic limit, the variance of the velocity field 
decays exponentially away from the wall, with a characteristic decay length given by the 
Stokes layer thickness = (21^/^1)^/^, where v is the kinematic viscosity of the fiuid, 
and f2 is the angular frequency of vibration of the boundary. Since the equations gov- 
erning the fiow are linear, we are able to obtain an analytic solution describing transient 
layer formation in the stochastic case, but only in the neighborhood of the white and 
monochromatic noise limits. We then show that for any finite correlation time the sta- 
tionary variance of the tangential velocity asymptotically decays as the inverse squared 
distance from the wall, in contrast with the exponential decay in the deterministic case. 
This asymptotic behavior originates from the low frequency range of the power spectrum 
of the boundary velocity. The crossover from power law to exponential decay is explicitly 
computed by introducing a low frequency cut-off in the power spectrum. 

We next investigate two additional geometries in which the equations governing fiuid 
flow are not linear, and show that several of the generic features obtained for the case of 



a planar boundary still hold. In the first case, we generalize the analysis of Schlichting 



(1979) concerning secondary steady streaming. He found that the oscillatory motion of 
the boundary induces a steady secondary flow outside of the viscous boundary layer even 
when the velocity of the boundary averages to zero. If the thickness of the Stokes layer, 
(5s, and the amplitude of oscillation, a, are small compared with a characteristic length 
scale of the boundary L {Sg L, a <^ L), then the generation of secondary steady 
streaming may be described as follows. Vibration of the rigid boundary gives rise to an 
oscillatory and nonuniform motion of the fiuid. The flow is potential in the bulk, and 
rotational in the boundary layer because of no-slip conditions on the boundary. The bulk 
flow applies pressure at the outer edge of boundary layer, which does not vary across 
the layer. The non uniformity of the flow leads to vorticity convection in the boundary 
layer through nonlinear terms. Both convection and the applied pressure drive vorticity 
diffusion, and thus induce secondary steady motion which does not vanish outside of the 
boundary layer. In the simplest case in which the far field velocity is a standing wave 
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U{x,t) — U{x)cos{ilt), the tangential component of the secondary steady velocity is, 



(1.1) 



where a; is a curvilinear coordinate along the boundary. In fact, Eq. (1.1) serves as the 
boundary condition for the stationary part of the flow in the bulk. Similar conclusions 
have been later reached by Batchelor (1967)| who studied sinusoidal oscillations of nonuni- 
form phase, and by Gershuni & Lyubimov (1998) who studied monochromatic oscillations 
of a general form. 



The second geometry that we address is the so called wavy wall ( Lyne (1971) ). The 
deterministic limit in which a wavy boundary is being periodically vibrated has been 
studied by a number of authors, mainly t o address the i nteraction between the flow above 
the sea bed and ripple patterns on it ( Lyne (1971)| , [Kancko fc Honjii (1979) , Vittori 
(1989) , Blondeaux fc Vittori (1994) and references therein). Lyne (1971) deduced the 
existence of steady streaming in the limit in which the amplitude of the wall deviation 
from planarity is small compared with the thickness of the Stokes layer. He introduced 
a conformal transformation and obtained an explicit solution in the limit of small kRe, 
where k is the wavenumber of the wall profile scaled by the thickness of the Stokes layer, 
and Re is the Reynolds number. The detailed structure of the secondary flow depends on 
the ratio between the wavelength of the boundary profile and the thickness of the Stokes 
layer. 

In Sections || and ^, we discuss how the results for these two geometries generalize to 
the case of stochastic vibration. Section || addresses the flow created by a gently curved 
solid boundary that is being vibrated randomly. The perturbation parameter that we use 
is the ratio between the amplitude of vibration and the characteristic inverse curvature 
of the wall. The ensemble average of the stream function is not zero, and hence there 
exists stationary streaming in the stochastic case as well. The average velocity diverges 
logarithmically away from the boundary because of the low frequency range of the power 
spectrum. We again introduce a low frequency cut-off tOc in the spectrum, and study 
the dependence of the stationary streaming on the cut-off frequency. We compute the 
stationary tangential velocity as a function of <C 1 and arbitrary /3, and find a weak 
(logarithmic) singularity as Uc 0. 

Section |^ discusses the formation of a boundary layer around a wavy boundary that is 
vibrated randomly. Positive and negative vorticity production in adjacent regions of the 
boundary introduces a natural decay length in the solution, thus leading to exponential 
decay of the fiow away from the boundary, even in the absence of a low frequency cut-off 
in the power spectrum of the boundary velocity. Steady streaming is found at second 
order comprising two or four recirculating cells per period of the boundary profile. The 
number of cells depends on the scaled correlation time Qt. 



2. Randomly vibrating planar boundary 

We first examine the case of a planar boundary that is being vibrated along its own 
plane. In this case the governing equations are considerably simpler then in the more 
general geometries discussed in Sections ^ and ^. In particular, the Navier-Stokes equation 
is linear, fact that allows a complete solution of the flow. Nevertheless, this simple solution 
still exhibits several of the qualitative features that are present in the case of random 
forcing by a curved boundary, namely asymptotic power law decay of the velocity field 
away from the boundary, and sensitive dependence on the low frequency range of the 
power spectrum of the boundary velocity. 
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Consider an infinite solid boundary located at z = 0, and an incompressible fluid that 
occupies the region z > 0. The Navier-Stokes equation, and boundary conditions are, 



dfU — vd^u, 



(2.1) 



u{Q,t) = uo{t), it((X),t) < cx), (2.2) 

where z is the coordinate normal to the boundary, u{z,t) is the x component of the 
velocity, and uo(t) is the prescribed velocity of the boundary. The solution for harmonic 
vibration uo{t) = Mocos(ni) was given by [Stokes (1851)| . It is a transversal wave that 
propagates into the bulk fluid with an exponentially decaying amplitude, 

u{z, t) = uoe~^/^'cos{VLt 

where 5s — (2j//ri)^/^ is the Stokes layer thickness. 



(2.3) 



2.1. Narrow band noise 



As discussed in the introduction, the main topic of this paper is to examine how the 
nature of the bulk flow changes when the boundary velocity uo(i) is a random process. 
Specifically, we consider a Gaussian process defined by 



K(t)) = 0, {u^{t)uM - W) e-l*-*'l/"cosl](t - t'). 



(2.4) 



This process is known as narrow band noise ( ^tratonovich (1967) ) . It is defined by three 
independent parameters: its variance (wq^, its dominant angular frequency O, and the 
correlation time r. Each realization of this random process can be viewed as a sequence 
periodic functions of frequency O, with amplitude and phase that remain constant for a 
time interval r on average. White noise is recovered when Jlr — > while D = {uq^') t 
remains finite, whereas the monochromatic noise limit corresponds to Qt — > oo, with 
(wo^) finite. Monochromatic noise is akin to a single frequency periodic signal of the same 
frequency, but with randomly drawn amplitude and phase. The relationship between the 
two can be illustrated by considering a deterministic function x{t) = xocos(flt) and 
defining the temporal average as, 



{x{t)x{t')) 



lim — 



dtx(t)x{t') = -^cos{n{t-t')). 



This average coincides with the ensemble average of the noise when (wq) 
power spectrum corresponding to the correlation function (2.4) is 



P{co) = 



2n 



(2.5) 
xl/2. The 



(2.6) 



dt UQ{t)e 



-iujt 



(2.7) 



We will also use the spectral density of the process uo{t), 

so that its ensemble average and correlation function are respectively given by 

(uo(w))=0, (Mo(t^)wS('^')) ='5(w-w')P(t^). (2.8a, &) 

We will often use dimensionless variables in which (uq) /il is the scale of P(w), and is 
the angular frequency scale. In dimensionless form. 



2n 



l + /32(tJ-l)2 l + /32(a; + l)2 



(2.9) 
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where /3 = fir. We have /.^ dwP(a;, (3) = 1, independent of /3, and also hm^^oo ^'('-<-') = 
— 1) + 6{uj + 1)] /2. Note that the power spectrum does not vanish at small fre- 
quencies. Instead, P(0, /3) = /3/(7t(l + /3^), which for large and small /? behaves as 
P{0,(3) ~ l/(7T/3) and P{0,(3) ~ /J/tt respectively. We will discuss separately the ef- 
fect of this low frequency contribution on the results presented in the remainder of the 
paper. 

2.2. Transient layer formation 

In the two limiting cases of white and monochromatic noise, it is possible to find an 
analytic solution for the transient flow starting from an initially quiescent fluid. The 
solution can be found, for example, by introducing the retarded, in finit e space Green's 
function corresponding to equation (2J_), with boundary conditions ( |2.2| ), 

1 



G{z,t\z\t') = 



t>t', (2.10) 



{iniy{t - t'))^^ 

and G{z, t\z', t') = for i < t'. If the fluid is initially quiescent, u(z, 0) = 0, we find 



with 



u{z,t) = i^ I dt'uo{t'){d,,G)^,^^ 



7? /4:U(t-t') 



(2.11) 



(2.12) 



[47Tl/3(i _i/)3]l/2 

Equations ( ^.ll| - pTT2 ) determine the transient behavior for any given UQ{t). 

If uo{t) is a Gaussian, white noise process, the ensemble average of Eq.(2.11) yields 
{u{z,t)) = 0. The corresponding equation for the variance reads 



{u^{z,t))^2D,.^ fdt' [{d.>G)^,J' 
Jq 



2Dv 



nz^ 



1 



z 

2iA 



(2.13) 



The variance of the induced fluid velocity propagates into the fluid diffusively. Saturation 
occurs for t 3> z"^ /2v, at which point the variance does not decay exponentially far away 
from the wall, but rather as a power law. 

{u\z,<x)) = ^. 2.14 

Ascertaining whether random vibration can induce flows far away from the boundary in 
more general geometries is one of the main motivations for this paper. 

Consider now the opposite limit of monochromatic noise with correlation function 



(uo(t)Mo(i')) = ("§)cos [n(t-t')]. 
Now using Eqs. ( ^J% and ( pl5| ) we find, ( |Carslaw fc Jaeger (1959)| ), 



2{ul) 



dcre 



d/i e 



261 VAi^ 



(2.15) 



(2.16) 



with K = z I (Avt)^!'^ . A closed form solution can only be obtained for long times. We 
find, 



2(ug) 



^/*^^sin(f^<-z/(5,) + C'(K^(t)). 



(2.17) 



2 711/2^2 

At long times, the variance propagates into the bulk with phase velocity {2v^)^l'^ , while 
its amplitude decays exponentially in space over the scale of the Stokes layer, and as 
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^-3/2 jj-^ time. In summary, the flow created by the vibration of the boundary propa- 
gates diffusively for white noise (z^ oc 2vt), and as a power law (z^ cx nvil^t^) in the 
monochromatic limit. 

2.3. Stationary variance for narrow band noise 

We have been unable to obtain a closed analytic solution for the transient evolution of 
the variance {u{z, i)^) when the vibration of the boundary is given by a general narrow 
band process. It is possible, however, to obtain the stationary variance of the velocity. 
Equation (|2.1[) can be rewritten in Fourier space as 



iuju{z,uj) = i'd'^u{z , uj) 



with 



u{z, t) 



Abj u{z, w)e" 



The bo unda ry conditions arc, u(0,lu) ~ uo(w), and u{z,uj) < oo at z 
of Eq. (|2.18| ) with these boundary conditions is 



u{z, LO) 



(l + isign(a;))(tj/2zy)i/2 



(2.18) 

(2.19) 
oo. The solution 

(2.20) 



We next choose as the time scale, and 6s as the length scale, and after some 
straightforward algebra we find. 



Iiz,P) 



POO 

/ dwP(w,/3)e 



(2.21) 



We have also used the fact that the power spectrum ( p. 91) is even in frequency. 

We next analyze the asymptotic dependence of /(z, 0) at large z. In this limit, /(z, 0) 
mainly depends on the low frequency region of the power spectrum; higher frequen cies 
are suppressed by the exponential factor. By using Watson's lemma ( Nayfeh (1981) ), we 
find, 

2/3 1 240/33(3/32 - l) i 



/(z,/3) 



0(z 



-10\ 



(2.22) 



7t(l+/32)z2 7T(1+/32)3 z6 

This asymptotic form at large z is valid for all /3. In particular, the dominant behavior 
for small and large (3 is 

2/3 1 _ 2 1 



I{z,f3) 



n z^ 



/(z,/3) 



7t/3z2 



respectively. Hence we recover the power law decay of Eq. ( 2.14 ) . 

We can also find the asymptotic behavior at large /3 that is uniformly valid in z, 

/(z, /3) = - ^ (^Ci(z) sin(z) - Si(z) cos(z) - i(e--Ei(z) - e^Ei(-z))^ + 0{f3-'), 

(2.23) 

where Ci and Si denote the integral sine and cosine functions, and Ei stands for the 



exponential integral function (Gradshteyn & Ryzhik (1980)). For z < 1, the variance 



decreases exponentially. At larger z, the exponential terms in (2.23) become small, so 
that the remaining asymptotic dependence for large z is given by Eq (p.22| ). The quantity 
J(z,/3) 02 computed both from (2.21) and the uniform expansion (2.23) is presented in 
Fig.|l|. For fixed (3, I{z, (3) asymptotes to 2/3/7t(l + 0^) outside of Stokes layer. This value 
is the coefficient of the leading term in the asymptotic expansion ( 2.22] ) . The expansion 
( p3| ) is a good approximation even for moderate values of /3. 



8 



D. Volfson and J. Vinals 



To summarize, the variance of the velocity field does not decay exponentially away 
from the wall for finite /3, but rather as the inverse squared distance. The crossover 
length separating exponential and power law decay increases with increasing (3. 

2.4. Low frequency cut-off in the power spectrum 



The coefficient of the leading term in (2.22) is in fact twice the value of P{0,P) 



/ 3/7t (1 + P'^). The algebraic decay of (wQ(z,t)) follows from the diffusive nature of Eq. 
(^.iD , and a non vanishing value of P{uj, (3) as, 13 ^ 0. Before we analyze in Sections ||and 
W how this behavior is modified by nonlinearities in the governing equations, we explicitly 
address here the consequences of a low frequency cut-off in the power spectrum. Of course, 
there always exists in practice a low frequency cut-off because of limited observation time. 
Furthermore, the low frequency range of the power spectrum of the residual acceleration 
field in microgravity {Vl/2n< 10^'^Hz) is fairly difficult to measure reliably. We therefore 
introduce an effective cut-off frequency in the power spectrum, ^ 1, and study the 
dependence of {u'^{z,t)) on lOc- The stationary value of variance of the velocity is now 
given by, 

^"'^:;y^ -/c(z,Ac..) = r ^.Pi.^P)e-"\ (2.24) 
By using Watson's lemma, we find for large z. 



I,{z,l3,L0c) = ^T-—^\^ + ^+\y.oX.\, (2.25) 




where h.o.t. stands for terms which are of higher order than terms retained under the 
assumption that both 1/z and Wc^^ are small but independent. For z ^ 1, but zujy^ <C 1 
the dominant term in ( |2.25 ) is 



^c(^,/?,^c)-^^^^^^, zul'^«l. (2.26) 

1 /2 • • • 1 /2 

On the other hand, if zujc ^ 1, the leading order term is now a function of — zujc 



Equations ( 2.26 ) and ( 2.27 ) show that at distances that are large compared with the 



thickness of the Stokes layer based on the dominant frequency 17, (v?{z^t)) decays al- 

1/2 

gebraically with z. There exists, however, a length scale zuic beyond which the decay 
is exponential. This new characteristic length scale is the thickness of the Stokes layer 
based on the cut-off frequency. This conclusion appears natural given the principle of 



superposition for the linear differential equation (2.1) 



3. Streaming due to random vibration 

Next we investigate to what extent the results of Section ^ hold in configurations in 
which the governing equations are not linear. We examine in this section the flow induced 
by a gently curved solid boundary that is being randomly vibrated. The boundary velocity 
is assumed to be described by a narrow band stochastic process, and hence our results 
will reduce to Schlichting's in the limit of infinite correlation time. However, for finite 
values of (3 the results are qualitatively different. The mechanism of secondary steady 
streaming generation described in the introduction is no longer valid because there is no 
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boundary layer solution at zcroth order. Vorticity produced at the vibrating boundary 
penetrates into the bulk fluid, at least perturbatively for small curvature. This results in 
a logarithmic divergence of the ensemble average of the first order velocity with distance 
away from the wall. A cut-off analysis is also presented, and similarly to that of Section 



2.4, it reveals the existence of an effective boundary layer of thickness based on the cut-off 



frequency. We then find an expression analogous to Eq. (1.1) as a function of /3 and uJc 
Define the following dimensionless quantities, 



{2{ul)y/^/nL, Rep ^ 2 (ul) /ili^, A ^ d? + / RepOj 



(3.1) 



Assume now that the characteristic scale of the boundary L is large so that e is a small 
quantity. If the Reynolds number Rep is assumed to remain finite, both conditions imply 
u/nL^ <g; 1. We next write the governing equations and boundary conditions in the 
frame of reference co-moving with the solid boundary and obtain for a two dimensional 
geometry (tildes are omitted), 

5.A^ + e^g^ = AV (3.2) 

o[z, X) 

= 0, d^^/j = at y = 0, (3.3a, b) 

d^tP = 2"^/^J7(x)mo(0 at z = oo, (3.4) 
where x, z are the tangential and normal coordinates along the boundary, and -0 is the 
stream function u = {dzip, —dxip)- We have also used the notation d{a,b)/d{z,x) = 
{dza){dxb) — {dxa){dzb) for the nonlinear term. The far field boundary condition is a 

1/2 

nonuniform and random velocity field, of the order of (uq) i with a spatially nonuniform 
amplitude U{x), and a stochastic modulation uo{t) which is a Gaussian stochastic process 
with zero mean, and narrow band power spectrum. We first expand the stream function 
as a power series oi e, ip ~ + etpi + . . . and solve ( |3.2| ) order by order. At order e° we 
obtain the following equation, 

{dtdl - dt)^o = (3.5) 

with boundary conditions, 

^po = 0, dzi' = at z = 0, (3.6a, b) 

dz^ = 2'^^^U{x)uo{t) at z = oo. (3.7) 

At this order, the equations effectively describe the fiow induced above a planar boundary 
with a far field velocity boundary condition that is not uniform. The solution can be found 
by Fourier transformation. We define, 

/■oo 

Mx,z,t)= duj Mx,z,uj)e'^'. (3.8) 



The transformed Eq. (3^) and the transformed boundary conditions (3^) allow separa- 



tion of variables. We define, 

Tpoix, z, uj) = 2"^/^C/(a;)uo(a;)Co(2;, w). 



so that Eq. ( |3.5|) leads to, 

{iujdl - d^Xo = 0, (3.9) 
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with boundary conditions = 0, dzC,Q = 0, at 2; = and dzQo = 1 at z 
solution is, 

= (l + isign(w))(|tj|/2)i/2). 



(0(2:, w) = -l/a + 2 + 1/ae ' 
At order e we find, 

{dtdl - dt)^. 



The 



(3.10) 



(3.11) 



with boundary conditions ipi = 0, 9z'0i = at ?/ = 0. The remaining boundary condition 
for ipi needs to be discussed separately. Consider first the classical deterministic limit 
which can be formally obtained by setting f3 — 00. Then, the right hand side of Eq. 
( 3.11 ) involves stationary terms (of zero frequency) , and sinusoidal terms with twice the 
frequency of the far field flow. Since the equation is linear, the solution ipi has exactly 
the same temporal behavior. In this case, it is known that it is not possible to flnd a 
solution for dzipi that vanishes at large z, but only one that simply remains bounded as 
z — > c«. By analogy, we introduce a similar requirement in the stochastic case of /? < 00. 
Since the zeroth order stream function diverges linearly, this condition simply amounts 
to requiring that the expansion in powers of e remains consistent. 

We now take the ensemble average of Eq. (3.11), and consider the long time stationary 
limit of the average {ip[^^ — limf^oo ii'i)), to find, 

dt4"^ = (a,Voa.a,Vo - d^^odl^o) ■ (3.12) 
The right hand side of this equation can be integrated from infinity to z. We obtain, 

ds)_UdU 



2 da; 



F(z,/3), 



(3.13) 



where 



and 



F{z,P) 



dujP{uj,/3)Q{cu,z), 



Qiio, z) = (-2 + 2Coa,Co - Coaf Co - Co^f Co). 



The power spectrum P{Lu,f3) is defined in Eq. (2_^). The constant that appears in the 
expression for Q{uj, z) comes from the pressure gradient imposed at infinity. 

We now proceed to solve Eq.( 3.13 ) subject to the boundary conditions V-'i"'' = c^zV'i*'' = 

(s) 

on the solid boundary, and dz'tpi < 00 as z — > 00. This is a boundary value problem 

(s) 

for ipi . In the limit /3 00, it can be solved analytically, and the result obtained by 
Schlichting is recovered, namely that the solution may be bounded at infinity simply by 
setting the homogeneous part of the solution equal to zero to satisfy the principle of 
minimal singularity (Van Dyke (1964)). Otherwise dzipi'^ is singular a fortiori. Since we 
cannot find an complete analytic solution for finite /?, we proceed as follows. We recast the 
boundary value problem as an initial value problem, and search for a boundary condition 
on d'^ip[^^ at z = so that the homogeneous part of the solution remains bounded. This 
boundary condition can be found analytically by integrating ( |3.13 ) from to z. We find, 

dLuP{uj,P)Q{uj,z'), 



dl4^\z)-dli,['\Q)^'^'^j\z' 



or after changing the order of integration. 



dli^['\z)-dl^[^\Q) 



U_dU_ 
2" d^ 



dujP{uj,P) 



dz'Q{uj, z') 



dz'Q{u 
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The second integral within brackets equals (l/2a;)^/^. In order to avoid a linear divergence 
of dzil^'^i^ {z) we equate the constant terms on both sides, thus obtaining the third initial 
condition 

= t^^^^((2(9 - P))'^' + (9 + 1)^/^ + (9 - 1)^/^), (3.14) 

where q = (1 + and the dependence of initial condition on (3 is shown explicitly. 

Equation (3.13) with its original boundary conditions, supplemented with Eq. ( 0.14| ) is 
an initial value problem, which we have solved numerically. 

Before presenting the numerical results, we study the asymptotic behavior of the so- 
lution for large z which is determined by the asymptotic form of F{z,P) at large z. By 
explicit substitution of the zcroth order solution we find 

Q(z, uj) = (-4cos(Z) - 2Z{cos{Z) + sm{Z)) + 2sin(Z))c"^ + 20^^^, (3.15) 

where Z — z{uj/2y/^. The leading contribution to F{z,0) as z — * cxo originates from the 
zero frequency limit of P(w,/3). Thus F{z,(3) ~ P(0,/3) dujQ(uj, z). The remaining 
integral may be easily calculated to yield the asymptotic form, 

6/3^ 1 

7t(l + (V) 

Therefore the stationary mean first order velocity m^*^ = dzip'i^ has a logarithmic asymp- 
totic form (see Eq. ( [3.13| )). This is to be contrasted with the power law decay of the 
stationary variance for the case of a randomly vibrating planar boundary. 

Finally, we have numerically calculated cJ^V'i'*^ (-2) for a range of values of (3. The results 



are presented in Fig. ^j. Equation (3.13) was integrated numerically with no-slip boundary 
conditions for ipi'\ and Eq. (3.14). The nmnerical results support our conclusion about 



the logarithmic divergence of u^^'^ with distance for any finite (3. As [3 increases the 
stationary mean velocity profile approaches a limiting form that corresponding to the 
monochromatic limit oi j3 ^ 00. In this limit we recover the Schlichting result, according 
to which dz'^p'^i\z) asymptotes to a constant value over a distance of the order of the 
deterministic Stokes layer. 

3.1. Low frequency cut-off in the power spectrum 

In Section ||, we showed that a low frequency cut-off in P{uj,(3) led to an exponential 
decay of the velocity outside an effective boundary layer of thickness determined by the 
cut-off frequency. We therefore examine here the consequences of a low frequency cut-off 
on the divergent behavior of the stationary average of the first order stream function. In 
order to find the asymptotic form of dzip[^\ we first integrate ( 3.13| ) twice from z = to 



z. By using the low frequency cut-off defined in Section (2^), we write 



dz4'\z,p,Lo,) = [Gciz,P,iOc) - G,(0,/3,c^,)] , (3.17) 

with, 

G,{z,l3,iu,)^ fdz' f dz'F,{z',P,u;,), 
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and, 

F,iz,P,u;,) = I dwP(c^,/3)Q(w,z). 



If we set Uc — and consider the monochromatic hmit of /? — ^ oo, we find that Gc(z, oo, 0) 
contains an exponential factor that vanishes at z ~ and that 00(0,00,0) = 3/4. 

Therefore the Schhchting result is recovered. An explicit form of Gc{z, f3, ojc) for any f3 and 
iVc can be obtained analytically, but it is far too complicated and we do not quote it here. 
It has a similar functional dependence as in the cut-off analysis of the planar boundary, 
and contains exponential terms involving (—20;^^). We find that dzip\^\oo, (3,uJc) = 
—UdU/2dx Gc{0, f3,ujc)- For finite but small ujc, we find, 

dA'\^:P:^c) = (2/3arctan(/3) - In ^^^^ + 0{u.l) (3.18) 

This asymptotic formula is compared with the numerically computed value of the tan- 
gential mean stationary velocity at large distances in Figure ^ Computations were done 
as described in the previous section, and for different values of ^ 1 and /3. In all 
cases dzip['^^ reached constant values at large enough z (the numerical value of infinity, 
Zoo, was chosen so that the change in velocity for z ^ Zoo was less than a prescribed 
tolerance). We also checked that any change in the boundary condition Eg. ( 3. 14 ) leads 
to a linear divergence in the tangential mean stationary velocity, thus confirming the 
adequacy of this boundary condition. The figure also shows that the computed values of 
dz'4'i'\z — 00, P,u!c) for small lOc are in a good agreement with (3.18). 

In summary, Eq. ( 3.18| ) shows that the tangential velocity away from the boundary 
asymptotes to a constant that is a function of /3, and has a weak (logarithmic) dependence 
on the cut-off frequency lUc- Therefore the asymptotic dependence in the stochastic case 
(with a cut-off) and the deterministic case is qualitatively similar, although the value 
of the asymptotic velocity of the former depends on (3. Note also that this asymptotic 
behavior only sets in for distances larger than (jz/wcO)^/^, a value that can be quite large 
in practical microgravity conditions. 



4. Randomly vibrating wavy boundary 

In the two previous sections we considered cases in which the characteristic longitu- 
dinal length scale of the solid boundary was either infinite or large compared with the 
characteristic amplitude of boundary vibrations. We examine here the case of a wavy 
boundary and study how comparable length scales in both directions infiuence the flow 
away from the boundary. In contrast with the Schlichting problem, the external applied 
flow is now uniform or, alternatively, the length scale over which the flow is not uniform 
is much larger that the wavelength of the boundary. Thus one anticipates that the nor- 
mal component of the flow that appears is caused by the wall profile. This flow interacts 
through nonlinear terms with the externally forced fiow that is parallel to the average 
profile of the boundary and gives rise to stationary streaming. Even for stochastic vibra- 
tion we show that positive and negative vorticity production in adjacent regions of the 
boundary introduces a natural decay length in the zeroth order solution, thus leading 
to exponential decay of the flow away from the boundary, even in the absence of a low 
frequency cut-off in the power spectrum of the boundary velocity. 

Consider a rigid wavy wall being washed by a uniform oscillatory flow parallel to the 
wall wave vector, 

u(a;,z = 00) = (uo(t),0). (4.1) 



Flow induced by a randomly ... 13 

We now assume that uo{t) is a narrow band Gaussian process. Assume also that the 
amphtude of the boundary profile / is small compared with both the Stokes layer Sg and 
the wavelength L, with Ss/L finite. The following dimensionless quantities are introduced, 

e = l/{v/9.f'^, Re=[2{ul) l^u]''^ , k = 2n{v /nf' / L, A = a| + 9| 

(4.2) 

referred to the Cartesian coordinate system sketched in Fig. ||. The solid boundary is 
located at 

fi{x) — fjee'x.piikx) + c.c. (4.3) 

with constant complex amplitude rj so that I?)! = 1/2. The dimensionless (and two di- 
mensional) Navier-Stokes equation (tildes are omitted in what follows) reads, 

d,Ai. + Re^^^^A^^, (4.4) 
o(z, X) 

with no-slip conditions at the boundary, 

tp — 0, dztp — at y = r]{x), (4.5a, b) 

and the imposed uniform flow at infinity, 

d^^p = 0, 9^-0 = 2~^/2mo(0 at z = oo. (4.6a, b) 

These equations depend only on three dimensionless parameters: e, the ratio of the am- 
plitude of the wavy wall to the boundary layer width; Re, the Reynolds number (the 
square root of Rep used in Section ||); and k, the wavenumber of the wall proflle in units 
of boundary layer width. We assume e ^ 1 and expand the stream function in a power 
series of e, 

= iZ-o + eV'i + • • • (4.7) 
The boundary conditions are likewise expanded in power series of e, 

ipiz)U=^ = Mz)\z=o + e(V'i(2)U=o + vd,Mz)U=o) + ■■■ (4.8) 



We now solve (4^) order by order in e. 

At zeroth order the wall is effectively planar. We decompose the Fourier transform of 
the str eam f unction as, '0o(a;,z,ci;) = 2~^/'^uo{uj)^o{z,uj). The function <^q{z,Lij) is given 



in Eq. ( [3.10 ). At this order, the solution is identical to that found for a planar boundary. 



At first order we seek a solution of the form, 

-01 = i) exp{ikx)(j){z, t) -f c.c. (4.9) 
so that the amplitude 4>{z,t) satisfies the Orr-Sommerfeld equation, 

{dtV~V^)(t>^ikRe{dlipQ-dziljQV)(t>, V = dl - k^ (4.10) 
with non-homogeneous boundary conditions, 

= 0, dz(t) = -d^^po at z = 0, (4.11a, &) 



= 0, dz(j) = at z = oo, (4.12a, &) 

The linear operator in the left hand side of Eq. (4.10) contains a significant difference 
with respect to that of Eq. (3.11), the equation governing the first order stream function 
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for the case of a slightly curved boundary. Both equations describe vorticity diffusion, 
but the biharmonic equation (4.10) contains a cut-off through the parameter k. Is is 
precisely this term that will lead to an asymptotic exponential decay of the velocity field 
sufficiently far away from the boundary for any finite [3. The exponential decay at long 
distances arises from the screening introduced by the simultaneous positive and negative 
vorticity produced at the troughs and crests of the wavy wall. 

In order to obtain a solution of the Orr-Sommerfeld equation (Eq. (4.10)), we further 

1 /2 

expand the amplitude (/>(z,t) in power series of kRe = 2n{2{ul)) ' im. This is the 
ratio between the amplitude of oscillation of the flow at infinity and the wall wavelength. 
We write, 

(j) = (j)Q + \_kRe(j)i + ... (4.13) 

The funct ion (j)Q obeys the linearized Eq. ( 4.10 ) with boundary conditions as in Eqs. 
( [4.1l[ - ^J2 ) with (j) replaced by 4>a. The Fourier transform of 0o is given by, 

U^^oj) - (2)-i/2^o(a;)-^ (e-"^ - e-^^) , (4.14) 

with p = (a^ + fc^)^/^, and the principal branch of the square root is assumed {^{p\ > 0). 
Recall that a = (1 + isign([j))(a;/2z/)^/^. The field 4>o describes vorticity diffusion near 
the wavy wall caused by the uniform but oscillatory far field flow. Both the spatial and 
ensemble averages of this flow are zero. However, the flow non-uniformity at this order 
induces mean flow at the next order, as it is readily apparent from the equation for 0i, 

(dtV - V^)(j)i = {dl^o - d,^o'D)^o, (4.15) 

with boundary conditions = 0, S^f/ii = at z = 0, oo. The field describes vorticity 
diffusion forced by the nonlinear interaction between (f>Q and tpQ. As was the case in 
Sec. 11 we focus on the long time limit of the ensemble average of Eq. ( [1.15| ) 0^*^ = 
Hmt^oo ^ X + c-c, where x is given by, 

V\ = ~^G{z,(3), (4.16) 

with, 

poo 

G{z,(3)= / dujP{u;,l3)Q{z,u;,(]), 
Jo 

Q(z, uj, 13) = a{p + fc)(2e-("*+'')" - e-("'+'^')" - e'"'). (4.17) 

The corresponding boundary conditions are homogeneous, x = 0, d^x = at z = 0, oo. 
The solution is, 



X{z,l3) = / dLuP{Lu,l3)xiz,^), 
Jo 

X{z, uj) = Aie-P'' + A2e-("*+'''^ -t- A3e-("*+'=)^ + (Bi + zB2)e-'=^ (4.18) 

where the functions A^, Bi depend on frequency and wavenumber, 

Ai^D/a^, A2 = D/{2a^p^), A3 = -D/{a^{a* + 2k)^ , D = a{p + k)/2, 

Bi = -(Ai +A2+ A3), B2 = {p- fc)Ai + {p + a* - k)A2 + a* A3. (4.19) 

Therefore the stationary part of the averaged first order stream function is given by, 

V'i^'' = ifci?er}exp(iA:a;)(x + X*) + C-C- (4.20) 



Flow induced by a randomly ... 15 

This solution shows that ip[^'' has a phase advance of n/2 with respect to the wah profile, 
and hence the flow in the vicinity of the boundary is directed from trough to crest 
{d^{x + X*)\z=o > 0). By shifting the coordinate system along the x axis we can change 
the phase of the wall profile so as to make it a simple cosine function. We consider 
ri{x) — eco8{kx) in what follows. 



Following Lyne (1971) , we now proceed to study the limits of k large and small, while 
kRe ^ 1. For k 1 the wavelength of the boundary profile is much smaller that the 
thickness of the viscous layer. In this case, a boundary layer appears of characteristic 
thickness 1/k. Screening between regions producing positive and negative vorticity oc- 
curs over a distance much smaller that the Stokes thickness based on the frequency of 
oscillation. The net vorticity does not diffuse even to distances of order z ^ 0(1), hence 
giving rise to exponential decay with an exp(— fcz) factor. The region in which the stream 

(s) 

function is not exponentially small depends on Z = kz. The explicit form of ipi may 



be obtained by direct expansion of the solution ( 4.20| ) in power series of 1/k, keeping Z 



fixed. The leading contribution to the steady part of the tangential component of the 
velocity is given by. 



■'"^ = kdz^j["^ - -^sm{kx) e-^Z{6 - Z^) (4.21) 



The boundary layer comprises two recirculating cells per wall period, located within 
< z < 1/k. The separation point is given by = 

In the opposite limit ol k <^ 1 one formally recovers the Schlichting problem in that 
the characteristic longitudinal length scale is much larger than the Stokes layer thickness. 
There is one fundamental difference, however, which can be seen from the solution, Eq. 



(4.18). It has two contributions: the first one is proportional Ai, arises from the particu- 



lar solution, and serves to balance the non-homogeneity in Eq. (4.16). This contribution 
decays within the Stokes layer. The second one is proportional to Bi, and arises from the 
general solution of the homogeneous part of the equation. This contribution decays over 
the stretched scale Z . It turns out that this second contribution introduces an additional 
separation of the composite boundary layer when < fc ^ 0.23. The lines of zero value of 
the longitudinal component of the steady velocity profile are shown Fig. |[ The location 
of the second separation point is entirely determined by that part of the solution that is 
proportional to Bi, and occurs at z ~ 0{l/k) » 1. 

The steady velocity may be obtained by expanding the exact solution (Eq. ( [4.18| )) in 
fii 

(outer solution u"i}). To leading order, we find. 



power series of k, first keeping z fixed (inner solution, u'l-^), and second keeping Z fixed 



u 



.;.,.,^-«...„,MU^(»»(='.-»»,.,,....,.,.leo,.,i, 



z' = z/2i/2 _ 0^^)^ (4 22) 



u['}{Z) -Rek'^sm{kx)[-{Z -l)e-^], Z - 0(1). (4.23) 



The solution for the inner and outer steady velocities was already obtained by Lyne 



(1971) by a conformal transformation technique. We further note that the inner and outer 
solutions can now be matched by requiring that u^^ (oo) = u [^J{0) = 3/ARek^sm{kx). 
Hence it is possible to construct a uniformly valid solution by adding the inner and outer 
solutions, and subtracting the first term of the inner expansion of the outer solution. We 
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find, 



4c (-z') Rek'^sin{kx) | [|-(sin(z') - cos(z')) + 2sin(z') + icos(z')] e^'^' 



+ iil + ^(2i/2^^'„ 1)6-2^/^/0.' I (4 24) 

We now turn to a numerical study of the case of finite [3. The boundary value problem 



(4.16) has been solved numerically by using a multiple shooting method for non stiff 



and linear boundary value problems (Mattheij & Staarink (1984)). The method has the 



advantage that the necessary intermediate shooting points are determined by the method 
itself, and that it can give the solution on a preset and nonuniform grid of points. The 



code was tested on the analytically known solution of the deterministic limit, Eqs. (4.18) 
and ( 4.19|). Our results are summarized in Figs. ^, ^ and ^ 



Figure g shows the separation points of the stationary velocity as a function of the 
boundary wavenumber k for a range of values of (3, including for reference the deter- 
ministic limit of /3 = oo (separation points are defined to be the zeros of the stationary 
tangential velocity). This figure shows that for fixed fc, the mean stationary velocity field 
may be comprised of two or four recirculating cells per wall period depending on /3. The 
first separation point is largely independent of [3, whereas the deviation of the second 
relative to its value in the monochromatic limit is proportional to /3/7t(l + /3^), the value 
ofP(0,/3). 

Our results in the limit fc <C 1 are presented in Fig. 0, where profiles of tangential 
component of the mean stationary velocity are plotted for fc = 0.1 and different values of 
(3. At large (3 (close to the monochromatic limit) the flow is comprised of four recirculating 
cells per boundary period. Upon decresing /3, the second separation point moves to infinity 
(see also Fig. ^), so that beyond some critical value of (3, only two recirculating cells 
remain. Further decrease in (3 results in the reappearance of the second separation point 
at infinity, which then continues to move towards decreasing z. The intensity of the 
recirculating modes does not change monotonically with (3 as we further discuss below. 
In the opposite limit of fc ^ 1 (Fig. ^ shows the case fc = 10.0; note that u"^^ is now 
normalized by Re/k^), the qualitative structure of the flow is largely independent of 13. 
The streaming flow has two recirculating cells per wall wavelength, and their intensity 
increases monotonically with decreasing /?. 

The complex dependence of the flow on (3 and fc can be qualitatively understood 
from the interplay between the width of the power spectrum (given by 1//3) , the viscous 
damping of each elementary excitation that depends on its frequency, and the penetra- 
tion depth of the flow field which is primarily dictated by the boundary wavelength. For 
small fc, large frequency modes are damped close to the boundary and do not penetrate 
much into the recirculating layers. Reducing /3 introduces high frequency components 
into the driving terms at first order, but they are dynamically damped. At the same 
time, the power in the dominant frequency components (around f2) decreases. Overall, a 
decrease in (3 then leads to a decrease in recirculation strength. As fc increases, larger fre- 
quencies contribute to the fiow over the entire range of the recirculating cells. Decreasing 
/3 decreases the strength of the dominant components, but increases the range of high 



frequencies that can contribute to the flow. From Eq. (4.16) one can show that the driv- 
ing contribution from higher frequencies which is contained in Q increases faster with 
frequency than the decreasing weight given to them by the power spectrum P{Lo,j3). 
Consequently, decreasing (3 (which amounts to moving towards the white noise limit) 
leads to increasing amplitude of the recirculation. 
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In summary, for any value of finite or infinite, the vorticity produced by vibration 
of the wavy boundary does not penetrate into the bulk farther than a distance of or- 
der of the wavelength of the boundary. However, there are qualitative differences with 
the deterministic limit in the character of the fiow within that layer. In particular, the 
structure and the intensity of the stationary secondary flow strongly depend on /3. 



5. Summary 

We have addressed the flow induced by a randomly vibrating solid boundary in an 
otherwise quiescent fluid. This analysis has been motivated by the random residual ac- 
celeration field in which microgravity experiments are conducted. The salient features of 
the flow are summarized below. 

When the solid boundary is planar, the flow field averages to zero (the average velocity 
of the boundary has been taken to be zero in all cases investigated), but its variance 
decays algebraically with distance away from the wall. This dependence follows from a 
non vanishing power spectrum of the boundary velocity at zero frequency. Introducing 
a low frequency cut-off in the power spectrum leads back to the classical exponential 
decay, with a rate that is determined by the cut-off frequency, Eq. ( 2.27| ). The amplitude 



of the decaying variance depends explicitly on the correlation time of the boundary 
velocity, (3 — fir, where is the dominant angular frequency of the power spectrum of 
the boundary velocity, and r is inverse spectral width (t is the correlation time of the 
boundary velocity). 

If the solid boundary is curved, steady streaming is generated in analogy with the 
classical analysis of Schlichting. The stationary part of the ensemble average of the sec- 
ondary velocity is nonzero, even though the boundary velocity averages to zero. In this 
case, we find that the leading contribution to the average stationary velocity diverges 
logarithmically with distance away from the boundary. In analogy to the planar case, 
the introduction of a low frequency cut-off in the power spectrum of the boundary ve- 
locity changes the asymptotic behavior quali tativ ely. The average stationary velocity 



asymptotes now to a constant, given by Eq. ( p.l8[ ). The asymptotic velocity explicitly 
depends on /3 and logarithmically on the cut-off frequency. This asymptotic behavior is 
not reached until a length scale of the order of the Stokes layer thickness that is based 
on the cut-off frequency. 

We have finally analyzed the case of a periodically modulated solid boundary in the 
limit in which the scale of the wall modulation is small compared to the thickness of the 
Stokes layer, and also when the spatial amplitude of the boundary oscillation is small 
compared with the wavelength of the wall profile. Cancellation of vorticity production 
over the wall boundary leads to exponential decay of the fluid velocity away from the 
boundary, with a decay length which is proportional to the wall wavelength, even if 
the zero frequency value of the power spectrum of the boundary velocity is nonzero. If 
the boundary wavelength is much larger than the Stokes layer thickness, we flnd steady 
streaming in secondary flow with two or four recirculating cells per wall period depend- 
ing on /3. On the other hand, if the wavelength is much smaller than the Stokes layer 
thickness, only two recirculating cells are formed regardless of the value of /3. Somewhat 
unexpectedly, the intensity of the recirculation can both increase or decrease with (3. 

This research has been supported by the Microgravity Science and Applications Divi- 
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Energy, contract No. DE-FC05-85ER25000. 



18 



D. Volfson and J. Vinals 



REFERENCES 

Alexander, J.I.D. 1990 Low-gravity experiment sensitivity to residual acceleration: a review. 
Microgravity sci. technol. 3, 52. 

Alexander, J.I.D. , Garandet, J. P., Favier, J.J. & Lizee, A. 1997 g-jitter effects on seg- 
regation during directional solidification of tin-bismuth in the mephisto furnace facility. J. 
Crystal Growth 178, 657. 

Alexander, J.I.D., Ouazzani, J. & Rosenberger, F. 1991 Analysis of the low gravity tol- 
erance of Bridgman-Stockbarger crystal growth. J. Crystal Crowth 113, 21. 

Batchelor, G.K. 1967 An Introduction to Fluid Dynamics. Cambridge: Cambridge University 
Press. 

Blondeaux, p. & ViTTORi, G. 1994 Wall imperfections as a triggering mechanism for Stokes 

layer transition. J. Fluid Mech. 264, 107. 
Carslaw, H.S. & Jaeger, J.C. 1959 Conduction of heat in solids. Oxford: Clarendon Press. 
DeLombard, R., McPherson, K., Moskowitz, M. & Hrovat, K. 1997 Comparison tools 

for assessing the microgravity environment of missions, carriers and conditions. Tech. Rep. 

TM 107446. NASA. 

Farooq, a. & HoMSY, G.M. 1994 Streaming flows due to g-jitter-induced natural convection. 

J. Fluid Mech. 271, 351. 
Gershuni, G.Z. & Lyubimov, D.V. 1998 Thermal Vibrational Convection. New York: John 

Wiley & Sons. 

Gershuni, G.Z. & Zhukiiovitskii, E.M. 1976 Convective Stability of Incompressible Fluids. 
Jerusalem: Keter. 

Gradshteyn, I.S. & Ryzhik, I.M. 1980 Tables of integrals, series and products. New York: 
Academic Press. 

Grassia, p. & HoMSY, G.M. 1998a Thermocapillary and buoyant flows with low frequency 

jitter. I. jitter confined to the plane. Phys. Fluids 10, 1273. 
Grassia, P. & Homsy, G.M. 19986 Thermocapillary and buoyant fiows with low frequency 

jitter. II. spanwise jitter. Phys. Fluids 10, 1291. 
Kamotani, Y., Prasad, A. & Ostracii, S. 1981 Thermal convection in an enclosure due to 

vibrations aboard a spacecraft. AIAA J. 19, 511. 
Kaneko, a. & HoNJii, H. 1979 Double structures of steady streaming in the oscillatory viscous 

flow over a wavy wall. J. Fluid Mech. 93, 727. 
Lyne, W.H. 1971 Unsteady viscous flow over a wavy wall. J. Fluid Mech. 50, 33. 
Mattiieij, R.M.M. & Staarink, G. W.M. 1984 An efficient algorithm for solving general linear 

two-point bvp. SIAM J. Sci. Stat. Comp. 5, 745. 
Nayfeh, A.H. 1981 Introduction to Perturbation techniques. New York: John Wiley & Sons. 
Nelson. E.S. 1991 An examination of anticipated g-jitter on space station and its effects on 

materials processes. Tech. Rep. TM 103775. NASA. 
Sciiliciiting, H. 1979 Boundary layer theory, 7th edn. New York: McGraw-HiU. 
Stokes, G.G. 1851 Trans. Camb. Phil. Soc. 9, 8, Mathematical and Physical Papers 3,1. 
Stratonovich, R.L. 1967 Topics in the Theory of Random Noise, vol. II. New York: Gordon 

and Breach. 

Thomson, J.R., Casademunt, J., Drolet, F. & Vinals, J. 1997 Coarsening of solid- liquid 

mixtures in a random acceleration field. Phys. Flmds 9, 1336. 
Van Dyke, M. 1964 Perturbation Methods in Fluid Mechanics. New York: Academic Press. 
ViTTORi, G. 1989 Non-linear viscous oscillatory flow over a small amplitude wavy wall J. Hydr. 

Res. 27, 267. 

Walter, H.U., ed. 1987 Fluid Sciences and Materials Sciences in Space. New York: Springer 
Verlag. 

Zhang, W., Casademunt, J. & Vinals, J. 1993 Study of the parametric oscillator driven by 
narrow band noise to model the response of a fluid surface to time-dependent accelerations. 
Phys. Fluids A 5, 3147. 



Flow induced by a randomly 



19 



0.15 - 




Figure 1 . Normalized variance of the t angen tial velocity for the case of a planar boundary com- 
pute d by numerical integration of Eq. (2.21) (symbols), and its uniform asymptotic expansion, 
Eq. ( ^.23| ), (solid lines). The function l(z, (3) asymptotes to a constant value outside of the 
classical Sokes layer based on fi. The uniform expansion remains a good approximation even for 
moderate /3. 
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Figure 2. Stationary first order velocity as a function of distance for a range of values of (3. All 
curves diverge logarithmically at large z, except for /3 = oo (monochromatic limit), in which the 
velocity asymptotes to a constant within the Stokes layer. This latter behavior reproduces the 
classical result of Schlichting. 




Figure 3. Stationary first order velocity as a function of distance for a range of values of (3. The 
power spectrum of the boundary velocity has a low frequency cut-off at ujc = 0.05. The velocity 
asymptotes to a constant that depends on the value of /3. 
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Figure 4. Asymptotic dependence of the stationar y vel ocity as a function of the cut-off fre- 
quency ujc- We show the case /? = 10 given by Eq. (3.18) along with the numerically obtained 
solution. 
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Figure 5. Schematic view of the geometry of the wavy wall studied in Section 




Figure 6. Separation points of the boundary layer over the wavy boundary as a function of its 
wavenumber k, for different values of The separation points are the loci of zero tangential 
velocity. At fixed k, the location of the second separation point depends strongly on /3. 
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Figure 7. Tangential component of the mean stationary velocity as a function of z for k = 0.1 
and a range of values of p. The case /3 = oo corresponds to analytic solution obtained by Lyne. 
The o ther curves are the numerical solutions of the boundary value problem defined by Eq. 
(4.16) and corresponding boundary conditions. 
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Figure 8. Tangential component of mean stationary velocity as a function of z for fc = 10 and 
a range of values of /3. The case (5 = oo corresponds to analytic solution obtained by Lyne. The 
other curves are the numerical solutions of the boundary value problem defined by Eq. (4.16) 
and corresponding boundary conditions. 



